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Abstract 

We study winding angles of oriented polymers with orientation-dependent 
interaction in two dimensions. Using exact analytical calculations, computer 
simulations, and phenomenological arguments, we succeed in hnding the vari- 
ance of the winding angle for most of the phase diagram. Our results suggest 
that the winding angle distribution is a universal quantity, and that the 9- 
point is the point where the three phase boundaries between the swollen, the 
normal collapsed, and the spiral collapsed phase meet. The transition between 
the normal collapsed phase and the spiral phase is shown to be continuous. 

I. INTRODUCTION AND SUMMARY 

The study of polymers is one of the most fascinating fields of current research, because 
of its relevance not only for material sciences, but also for the understanding of proteins. 
Depending on the chemical and physical environment, a polymer in a dilute solution can be 
either swollen or collapsed, or at the ^-point, which is the boundary point between the two. 
Such polymers can be modeled by interacting self-avoiding random walks (SAWs) with 
an interaction energy e between (non consecutive) bonds that are on the same plaquette 
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of the underlying lattice. As the temperature is decreased, the SAW undergoes the above- 
mentioned transition at a the 6'-temperature, provided that e < 0. The value of the exponent 
V that characterizes the relation between the polymer length (monomer number) and its 
radius of gyration R/vs, v = 3/4, 4/7, and 1/2 above, at, and below the ^-temperature in 
two dimensions. 

The phase diagram for the polymer collapse becomes more complex when the polymers 
are oriented^ i.e., when they look different in the two directions along the chain, as, e.g., 
for A-B polyester 0. In this situation, the interaction energy between nearby monomers 
depends in general on whether their relative orientation is parallel or antiparallel. When 
the attractive interaction between parallel monomers is sufficiently strong, the collapsed 
polymer winds up to form a spiral. A phase diagram, based on numerical work and exact 
results, was suggested in 0. It contains three phase boundaries, separating the swollen, the 
normal collapsed, and the spiral phase (the latter being also a collapsed phase), and meeting 
at one point. 

In that phase diagram, the line along which parallel and antiparallel interactions are 
equally strong plays no special role. If this is correct, the values of the critical exponents 
do not fall into universality classes that are determined by symmetries. In fact, conformal 
field theory suggests that the exponent associated with the partition function (usually 
denoted 7) may depend continuously on the parallel interaction energy in the swollen phase, 
while the exponent v remains constant. This supposed nonuniversality of 7 in the swollen 
phase is complemented by the observation that at the collapse transition 7 assumes different 
values on different lattices: Its value on the square lattice is different from that on the so- 
called Manhattan lattice, where each bond of the lattice has a preassigned orientation, 
thus naturally excluding parallel contacts. Numerical studies of the swollen phase 
however, show no indication of a variation of 7 with the parallel interaction strength. These 
studies involve exact enumeration Monte Carlo simulation 0,^, and transfer matrix 
calculation on a strip of finite width ||^. In all these studies, as well as in our own exact 
enumerations, the (very small) variation of 7 with the parallel interaction strength decreases 
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with increasing polymer length, making it unlikely that 7 should show nonuniversal behavior 
for much larger polymer length. 

Furthermore, transfer matrix calculations on the collapse line suggest that, whenever 
the antiparallel interactions energy is lower than the parallel one, 7 has the same value as 
on the Manhattan lattice, and that the ^-point (where parallel and antiparallel interactions 
are equally strong) is the point in the phase diagram where the three phase boundaries 
meet. If this scenario is correct, the exponent 7 is a universal quantity, and its value at the 
collapse transition depends only on whether the symmetry between parallel and antiparallel 
interactions is broken. 

In this paper, we study the winding angle distribution for polymers in two dimensions 
with orientation-dependent short-range interactions. Two monomers that have a parallel 
contact are connected by a loop that encloses one of the end points of the polymer, i.e. 
their winding angle differs by 2tt. Since the winding angle is so closely related to the 
occurrence of parallel contacts, it should be equally sensitive to a change in the parallel 
interaction strength as the exponent 7. In fact, we find analytically a different winding 
angle distribution at the collapse transition on the Manhattan lattice and on the square 
lattice. As for the exponent 7, the question arises whether the winding angle distribution 
is determined by simple universality criteria. Usually, winding angle distributions depend 
only on universal features like symmetries and interaction range, when the length of the 
polymer is sufficiently large 0. The important result of this paper is that the winding 
angle distributions for oriented interacting polymers are also universal. The main evidence 
comes from the collapsed phase. Minimizing the free energy, we find that the winding angle 
distribution in the collapsed phase depends only on whether the symmetry between parallel 
and antiparallel contact energies is broken. Just as the transfer matrix calculation of [Q, 
this modifies the phase diagram suggested in and moves the phase boundary between the 
normal collapsed and the spiral phase on the symmetry line, where parallel and antiparallel 
interactions are equally strong. Our argument shows also that the phase transition along 
this line is continuous, in contrast to |@J^. The 6'-point is the point where the three phase 
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boundaries meet. Due to its special role, it can have an exponent 7 and a winding angle 
distribution that is different from all other points, without leading to nonuniversal behavior. 

Certainly, this phase diagram is very appealing due to its simplicity. Remarkably, we 
arrived at our results independently of 0, and prior to learning about that work. While 
1^ provides the stronger numerical evidence for the special role of the 6'-point, our work 
gives insights in the underlying physics of oriented polymers: Whenever the antiparallel 
interaction energy is more negative than the parallel one, the winding angle of the polymer 
decreases during the collapse. It is confined in the collapsed phase and becomes zero in 
the ground state, where the end points must be at the surface. When the two interactions 
are equal, however, the variance of the winding angle increases during the collapse, and is 
always proportional to InA^. When the parallel interaction is stronger than the antiparallel 
one, the collapsed phase has an overall spiral shape, the number of parallel contacts being 
proportional to the length of the polymer. 

The outline of this paper is as follows: In section we describe some basic properties 
of self-avoiding interacting walks, introducing and discussing the partition function and the 
general form of the winding angle distribution. In section |T|, we derive an exact expression 
for the winding-angle distribution at the collapse transition in the absence of parallel con- 
tacts between monomers. We also show results of a Monte Carlo simulation that agree well 



with this analytical result. In section we conjecture the variance of the winding angle 
for most of the phase diagram, and we argue that the transition from the normal collapsed 
to the spiral state occurs when parallel and antiparallel interactions become equally strong. 
We compare to the phase diagram obtained from exact enumeration and explain the origin 
of the discrepancies. Section summarizes and discusses our results. 

II. GENERAL PROPERTIES OF INTERACTING SELF-AVOIDING WALKS 



4 



A. The partition function and the critical exponent 7 



We model a polymer in two dimensions by a SAW on a two-dimensional square lattice, 
as shown in Fig. |I]. (In some cases, also a hexagonal lattice is chosen.) The steps of the 
walk, which coincide with lattice bonds, can be viewed as monomers. Each site of the lattice 
can only be visited once. This condition models the excluded volume effect of the polymer. 
Starting at one end point and stepping along the trajectory of the SAW, each bond can 
be assigned the direction in which it is passed. By this procedure, the polymer obtains an 
orientation. Short-range interactions between monomers are taken into account by assigning 
an energy or to each pair of (non consecutive) monomers that lie on the same plaquette 
and that have an antiparallel or parallel relative orientation. The weight for each parallel 
(antiparallel) contact is given by Up = exp(— ep/fc^T) {ua = exp{—ea/kBT)), where T is the 
temperature and fc^ is the Boltzmann constant. 

The partition function for a polymer of N steps is then 

Zn= gN{ma,mp)u:^^u';'' , (1) 

ma,mp 

where ma and nip are the number of antiparallel and parallel contacts, and gNi'^ai^p) is 
the number of configurations (starting at a given point) with these contact numbers. 

For oja = ojp = one has a normal SAW without any interaction except the self- 
avoidance, and the partition function is identical to the total number of SAWs starting at a 
given point, which is known to be 

Zn ~ A/x^AT^-i , (2) 

with A 1.771, yU ^ 2.638 and 7 = 43/32 P|. This form of the partition function is 
believed to hold also for Ua^Up 7^ 1, as long as the polymer is in the swollen phase, with 
different values of A and /i. 

The mean number of antiparallel and parallel contacts in the swollen phase is 

(9 In Ztv 

{nii) = uJi— 

ouji 

5 



\ OUi OUli OOJi j 

with i = a or i = p. While A can depend on both Up and Ua, the free energy per monomer 
fi does not change with ujp. This was exactly proven in for tUp < Ua with c^a = 1, and 
— assuming the scaling form Eq. (0) — it implies that the number of parallel contacts 
increases not faster than logarithmically with A^. This is not surprising when one realizes 
that a SAW can have antiparallel contacts anywhere along its trajectory (i.e., {rria) oc A^), 
while the average number of parallel contacts should not increase faster than the number of 
windings, which in turn increases logarithmically in A^. 

On a more phenomenological level, one can make the following entropy consideration: 
as long as Up < Ua, there is no energetic disadvantage for the polymer to have its endpoints 
at the surface, in which case it has no parallel contacts. Such configurations correspond 
to the case ujp = 0. We can therefore assume that the probability that a randomly chosen 
configuration (for oOa > ujp > 0) has both end points at the surface is not smaller than the 
square of the ratio of the surface area to the volume, proportional to oc N^^". The 
entropy loss per monomer due to the restriction of the endpoints to the surface is of the order 
21/111 N/N and vanishes in the thermodynamic limit A^ oo, and so does the change in free 
energy per monomer. This argument is not restricted to the swollen phase. Consequently, 
the phase boundary between the swollen and the collapsed phase cannot depend on Up. 

When Up is increased beyond Ua, there must be a point where the free energy becomes 
dependent on Up, as proved in 0. This point, which is a nonanalyticity of the free energy, 
marks the phase transition to the spiral phase (if one assumes the simplest scenario of only 
one phase transition along a line Ua =const). The qualitative phase diagram is shown in 
Fig. |. 

So far, we have not yet discussed the last term in Eq. (|^). It seems implausible that 
a critical exponent, which is a universal quantity, should vary within one phase, and this 
term should therefore vanish. In particular, a repulsive interaction between monomers (i.e., 
Up^a < 1) is nothing else than an increased excluded volume, which can hardly modify the 
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value of a critical exponent. If d'^/duop vanishes, the mean number of parallel contacts is 
independent of for large A^, i.e., it saturates, in agreement with the above-mentionend 
numerical results [^1]-^. If d'y/dujp did not vanish, as suggested by the conformal field 
theory , the number of parallel contacts would increase logarithmically in A^. An increase 
in the excluded volume for parallel contacts, i.e., an increase in Up by 6ujp, would then have 
a similar effect as an increase in polymer length from A^ to A^ + a In A^ in the expression for 
Zn above, with a = {d'y/du!p){6u!p/\rLfi). 

The following argument suggests that the number of parallel contacts saturates in the 
limit A^ — > oo. We assume that the polymer is radially scale invariant, i.e. that it is 
statistically mapped onto a polymer of length bN, when the coordinates r of all monomers 
are scaled to b'^f (the starting point of the polymer being the origin of the coordinate system). 
This is equivalent to the statement that, when represented in the ^ — Inr-plane, the polymer 
is translation invariant in Inr-direction. For each value of r, there exists consequently the 
same mean number of parallel "close encounters" over a given distance A(lnr). However, 
only for A(lnr) = 1/r such a close encounter corresponds to a parallel contact. The total 
number of parallel contacts is then proportional to J^^^^ {l/r)d{\nr) and converges for 
A^ — > oo. This argument is in agreement with the numerical observation in ^ that the 
total number of "loops" saturates for large A^. We conclude that the universality hypothesis 
holds also for the exponent 7, and that 7 is constant in the swollen phase, provided that 
our assumption of radial scale invariance is correct. It might be that the conformal field 
theory in does not capture correctly the difference between parallel contacts and "close 
encounters" . 

At the 6'-point Up = Ua = 000, the value of 7 (denoted 74) is known to be 74 = 8/7 [1^. It 
is also known that 7^ = 6/7 for a self-avoiding walk that is part of a percolation cluster hull 



and therefore has no parallel contacts ||Tl|. Since v = A/7 JT^ for these walks, they are at 
the collapse transition. The authors of argue that 74 = 6/7 for all Up < Ua = uj0. If this is 
correct, the partition function contains a crossover term close to the ^-point Up = Ua = 000, 
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Zm ^ fi^'N'/'fiAiuN'') , (4) 

where /(x) ~ const, for small x, and f{x) oc x'^^"^^ for large x, and Au; = Ua — ojp. is a 
crossover exponent. A calculation analogous to Eq. (^) gives then {nip) oc at the ^^-point. 
Since the mean number of parallel contacts saturates in the swollen phase and increases as 
in the collapsed phase (see below), we expect < \1' < 1/2. For a related problem, 
the adsorption of a self-interacting polymer at a surface, the corresponding exponent has 
the value ^ = 8/21 ||12|. Since the adsorption of a polymer at its own surface (i.e., spiral 
formation), is somewhat different, the two crossover exponents need not be the same. 

In the low-temperature phase, the polymer has a finite density. Therefore u = 1/2, and 
surface effects become important. The partition function is assumed to have the general 



form ||T3| 



Zj^^Afi^K'^N^'-'. (5) 

We will argue below that for ujp = LUa the number of parallel contacts is proportional to V^, 
while it saturates for Up < Ua- 

B. Winding angle distributions 

For self-avoiding walks, the winding angle distribution is generally described by a Gaus- 
sian 

p{e) oc exp(-^V2C In A^) , (6) 

with a variance Cln |[l^,[15|. Such Gaussian distributions with a variance proportional to 
In A^ occur generically for radially scale invariant polymers, given that the winding center is 
visited only a finite number of times ||^. The value of C is C = 2 in the swollen phase and 
C = 24/7 at the ^-point of a polymer with no orientation dependence in the interaction ||15|| , 



as obtained from an analytical calculation on an hexagonal lattice. The value C = 2 is nicely 
confirmed by an exact enumeration on a square lattice, where we obtain C = 2.0005(6) for 
polymers up to length A^ = 26, using differential approximant analysis ||16|| . 
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The value of C in the collapsed phase (for ojp = Ua) is not known, but it is larger than the 
previous two values, since the winding angle apparently increases during the collapse. For 
dense SAWs, C is known to be C = 4 Since they have no self-interaction apart from 
self-avoidance, dense polymers have their finite density due to an external pressure and are 
"hot" , in contrast to collapsed polymers, where the attraction between monomers determines 
the density. For this reason, the density of monomer-monomer contacts is different in both 
cases, and it is not clear whether C can be the same. 

In the next section, we will derive C = 6/7 for a self-avoiding walk that is part of a 
percolation cluster hull and therefore has no parallel contacts. Since u = A/7 for these 
walks, they are at the collapse transition. The value of C in most other parts of the phase 



diagram will be discussed in Sec. IV 



III. THE WINDING ANGLE DISTRIBUTION AT THE COLLAPSE 

TRANSITION 

In this section, we calculate the winding-angle distribution for a SAW on a Manhattan 
lattice at the collapse point. Fig. ^ shows such a walk. The Manhattan lattice is an array of 
alternating one-way streets, thus not allowing parallel contacts and always keeping the way 
back to the origin open. SAWs on this lattice can be grown kinetically in a very efficient way, 
since they get trapped only through loop formation. One starts at the origin and constructs 
a path by going at each step in one of the allowed directions. When the path closes to a 
loop, it is cancelled. Since this procedure gives contacts a higher statistical weight than free 
steps, the polymer has an effective interaction. In fact, one can show that the path can be 
mapped onto the perimeter of a percolation cluster [|1^ , which in turn is known to have an 



exponent v = A/7 |T^. This means that the polymer is at the collapse point. 

In order to find its winding angle distribution, we have to calculate the winding angle 
distribution of parts of percolation cluster hulls. The procedure is similar to the one in 
r^,|T^, and is conveniently performed on an hexagonal lattice. We start with the 0(n)-loop 
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model |T8[ with the partition function 

Zo(n)= jX{d''SkX{{l + (3S,-S,) 

= Y.P'nr (7) 
g 

i,j,k are lattice sites, nearest neighbors, and S an n-vector: = n. In the second 
line, the sum is performed over all graphs Q of P non-intersecting polygons of total length 
/. For n G [—2, 2], the loop model has a critical point Pc = [2 + (2 — n)^/^]^^/^. For n = 1 
and P = I, each loop has the same weight. These loops can be interpreted as percolation 
cluster hulls for site percolation on the dual triangular lattice at the percolation threshold 
Pc = 1/2. We are interested in the winding-angle distribution of a segment of a loop. Fig. ^ 
illustrates the following calculation. We look at the function 

Go{n){X - Y, d, ea) = ^o(iV)(^i) exp[zei7r(ni + n[) + ie27r{n2 + n'^)] (8) 

gi 

Qi are like the graphs Q above, but with both points i and j lying on the same loop. Both 
parts S and S' of the loop are given the same orientation from i to j. rii (ri'^)and n2 (n'g) 
count the number of intersections of the oriented paths S (S") with Li and L2 respectively, 
crossing in different directions having a different sign. Without the phase factor, Go(n) is a 
four-spin correlation function, with two spins at i and two spins at j. 

To calculate Eq. (|^), we transform it into a solid-on-solid (SOS) model. Height variables 
are defined on the centers of the hexagons, such that two adjacent heights are equal or 
different by ivr. The polygons, once arbitrarily oriented, become domain walls with a step 
of +7r on the left of any oriented line. Along the straight line connecting X to Y , the height 
changes by 27r. In the SOS language, this corresponds to a dislocation line with a vortex 



at X and an antivortex at Y ||T^. At the vertices of the honeycomb lattice, the SOS walls 
turn by ±7r/3. The SOS weight Wsos is calculated as the product along the walls of local 
factors /3exp(m) (/5exp(— m)) at each left (right) turning vertex. Summing over the two 
independent orientations of each polygon (except, of course, the special polygon connecting 
i and j) gives a phase factor 2cos6'U for each polygon. The SOS weight of a graph Qi is 
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then 

W"sos(^i) = P^n^ exp[iu{n+ - n_) + iu{n'_^ - n'_) + 6m(P+ - P_)] . 

n+ and [n'j^ and n'_) are the total number of local left and right turns of path S (5"). 
P+ and P_ are the total number of right and left polygons surrounding i and j. In the 
asymptotic limit, we have n+ — ri_ = 6{ni + 722) and — n'_ = 6{n[ + n'g). Let us define 
the SOS correlator 

GsosiX - Y, e'l, e'2) = E W^sos(6?i) exp(2e;0(X) + ze',(P{Y)) . 

Equating this correlator with Eq. (Q) above gives n = 2 cos 6u and e'l = ei + eo, 63 = 62 + eo, 
ei + 62 = 0. The new constant eo is eo = —Gu/n. 

Now, at /3 = 1, the SOS model renormalizes onto the low-temperature phase of the 
Coulomb-gas model for g = 2/3 {ii n = 1) |T^, and 



Gsosle'i, e'2) = \X- y |44/9+9'^i™2 . (9) 



The magnetic charges mi and m2 are mi = —mi = 1/2 ||T^, due to the vortex pair. 
The winding angle is finally extracted from 

(exp[ie7r(ni - ^2 + n[ - n2)])o{N) = exp{-e'^g~^ In \X - y|) . 

Fourier transformation yields immediately a Gaussian distribution for rii — ^2 + n[ — n'2. 
Each of the two paths S and S' has the same number of intersections with L and L\ and 
therefore ni = n[ and n2 = n^. (In certain situations, ni and n'-^ differ by one, but this effect 
can be neglected in the thermodynamic limit, where the number of intersections becomes 
very large.) In terms of the angle 6 = 7r(ni — 77,2 + n[ — n'2) the resulting distribution reads 

P{e) = {IQng^^ In \X - Y\)-^^^ exp{-ge^/16 In |X - f |) . 

Since both paths make exactly the same contribution to 6, the winding angle of one path is 
given by 6/2. For large distances \X — Y\, the windings around X and Y are independent 
from each other and have the same probability distribution. Replacing |X — y| by N'^ [N 
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being the length of the polymer), we arrive at the winding angle distribution of path 5* 
around point X, 

P{e) = {Ani^g~HnN)-^/'^exp{-ge'^/2iy\nN) . (10) 

Inserting g = 2/3 and u = 4/7, we find C = 6/7. 

Our numerical simulations confirm nicely this result. Fig. |^ shows the winding angle 
distribution for polymers of length up to 10^, and the variance of the winding angle for 
length up to 10^. The solid line is the analytical result. 

In the swollen phase, we believe that the constant C is C = 2 for a SAW on the Manhattan 
lattice, i.e., the winding angle distribution is the same as for the normal SAW. The above 
calculation cannot be repeated in the swollen phase, since /5 < 1 in Eq. (|^) above. Although 
there exist still two paths connecting i and j, they now have different weight. A path of 

steps will close to a loop after a mean number of steps that diverges, and therefore the 
second path has in fact no weight at all. We expect to arrive at the same situation as for 
a normal SAW, which is discussed in [|T^]. Conformal field theory suggests also that C 
does not depend on the strength of the parallel interactions. Only the magnetic charges in 
Eq. are affected by it, and these drop out when the winding angle is calculated. 

The calculation of this section can easily be generalized to "watermelon configurations" , 
where the points i and j are connected by L paths. The constant C characterizing the 
winding angle distribution for one of these paths is then Cl = 'iv/ L'^g. 

IV. WINDING ANGLES FOR COLLAPSED POLYMERS, AND THE PHASE 

DIAGRAM 

The above result C = 6/7 was obtained in a situation where the polymer cannot make 
any parallel contacts. We therefore suggest that this corresponds to the case Up = 0, c^a = ue- 
However, a SAW on a square lattice that has no parallel contacts is different from a SAW 
on a Mahattan lattice, since it can get trapped without loop formation. Only if the range 
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of the repulsive parallel interaction is extended to next-nearest neighbors, the way back to 
the origin remains always open. If one assumes that all these situations are equivalent to 
each other, one must draw the conclusion that the precise range and form of the repulsive 
interaction between parallel bonds is not important, and that C = 6/7 holds on a finite 
part of the transition line. Assuming only one nonanalyticity on this line, we conclude that 
C = 6/7 on the entire line a;^ < a;a, and that it jumps to C = 24/7 at the 6'-point. The 
^-point is consequently the point where the three phase boundaries meet. 

Of course, other interpretations are in principle possible, e.g., that the interacting SAW 
on the Manhattan lattice falls into a separate universality class. Our argument, however, is 
supported by our exact enumeration data. Although they do not allow a good estimate of 
the value of C, they show clearly that for both cases cUp — 1 and cUp — the winding angle 
decreases during the collapse (i.e., with increasing cUa), while it increases during the collapse 
along the line cUp — cua- The constant C should therefore be smaller than 2 for ujp < uja — oJe-i 
and assuming universality, the conclusion C = 6/7 follows naturally. It is the equivalent of 
the transfer matrix result that the exponent 7 assumes its Manhattan lattice value 6/7 for 
all ujp < uJa — oj0, and that it jumps to 8/7 at the ^-point ujp — uJa — oJe- 

Our universality conjecture finds its strongest support when one studies the collapsed 
phase of the polymer. In the following, we argue that C = in the collapsed phase whenever 
ojp < oJa {C — Q should also hold for a collapsed polymer on the Manhattan lattice), and 
C — 00 when ujp > cUa- We start with the assumption that the winding angle distribution 
does not change when one goes slightly to the right or to the left of the line ujp = uJa, and 
we lead this assumption to a contradiction. 

On the diagonal cUp — cUa, the winding angle distribution is Gaussian with some unknown 
but finite constant C. Since C increases when going from the swollen phase to the ^-point, 
it certainly becomes even larger in the low-temperature phase. A polymer in the collapsed 
phase on the diagonal cUp — cUa has a finite density of contacts along its trajectory, leading 
to dfUp oc Rd6, and (with dO (x dlnN and R oc iV^/^) to cx iV^/^. 

Let us first consider the case cUp > cUa, with cUp — cUa — Ao; <^ 1. We compare the free 
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energy of a collapsed polymer in a globule configuration similar to the one on the line Up = Ua 
to the free energy of a spiral that is composed of globules of n monomers (see Fig. ^). The 
difference in internal energy between the two is (neglecting constant coefficients) 



AU ~ TAlu [VN - {N/n)y/^ 

When one transforms a globule to a spiral, one breaks 0{-\/N) parallel contacts, and one 
creates 0{N^/n) new parallel contacts. The leading contribution to the difference in entropy 
is 

AS ~ -{N/n) Inn. 

This is the number of globules times the entropy loss per globule when the end points of the 
polymer are restricted to the surface of the globule (see Sec. LI A ). The entropy difference 



between one large globule of monomers and N/n globules of n monomers (without any 
constraint for the end points) increases slower than 0{N), since the entropy is an extensive 
quantity. For any Auj, the gain in binding energy —U is larger than the loss in entropy, 
when n is sufficiently large. By minimizing the free energy AF = AU — TAS, we find (to 
leading order) 

nocAu;"^ (11) 

The correlation length ^ oc ^/n, which is proportional to the globule radius, diverges as 
^ oc (Au;)"^. These results are correct for N ^ n ^ 1. 

Spirals have a considerable entropy close to the transition. However, it is difficult to see 
the continuous character of this phase transition in simulations. The spiral shape can only 
be seen when polymer length is much larger than the globule size n. The transition appears 
to be shifted to the right by a distance Alu = 0(l/-\/iV), as in 0. 

Now, we consider the case Up < Ua, with Ua — Up = Alu <^ 1. We compare the free energy 
of a collapsed polymer in a configuration similar to the one on the line ujp = Ua to the free 
energy of a polymer that has its end points at the surface. Bringing the end points at the 
surface, replaces 0{^/N) parallel contacts by antiparallel contacts, decreasing the internal 
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energy by an amount proportional to AU ~ Auj\fN . The entropy loss due to this restriction 
is of the order In A^, as shown previously. Having the end points at the surface means always 
a decrease in free energy in the thermodynamic limit N ^ oo. This means that the winding 
angle and the number of parallel contacts must saturate in the thermodynamic limit. Only 
an initial and finite segment of length n oc (Au;)^^ of the polymer may behave like a polymer 
at the transition line. As for the spiral phase, the correlation length diverges as ^ oc {Au)^^. 

We thus have shown that the constant C of the winding angle distribution in the collapsed 
phase has the same value C = for all Up < Ua- It would be rather surprising if this universal 
feature did change at the collapse transition or in the swollen phase, where contacts play a 
less important role. The study of the collapsed phase therefore provides convincing support 
for the hypothesis that the winding angle distribution is a universal quantity. The constant 
C characterizing the winding angle distribution can now be given for most of the phase 
diagram, as indicated in Fig. |^. 

We have also seen that the number of parallel contacts saturates in the collapsed phase 
for Up < uja- This result should also hold in the swollen phase, where contacts are less 
important, leading to a universal exponent 7. 

V. SUMMARY AND DISCUSSION 

In this paper, we have studied oriented polymers with orientation-dependent interaction. 
We have argued that both the winding angle distribution and the exponent 7 are universal 
quantities, in agreement with transfer matrix calculations for the exponent 7 0. This result 
is closely tied to the observation that the 6'-point is a special point in the phase diagram, 
where three phase boundaries meet. When parallel contacts have a different energy than 
antiparallel contacts, a symmetry is broken, and a phase transition takes place. In contrast 
to 0, we predict that this phase transition between the collapsed and the spiral phase is 
continuous. The existence of this phase transition is particular to two dimensions, since in 
higher dimensions a parallel contact can locally be transformed into an antiparallel contact, 
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without changing the conformation of the polymer at a large scale. 

We succeeded in obtaining the winding angle distribution for the case where antiparallel 
contacts dominate, at the collapse transition as well as in the low-temperature phase. We 
have also argued that the number of parallel contacts saturates in the thermodynamic limit 
whenever the antiparallel energy is larger than the parallel one. 

There are three challenges left: The winding angle distribution along the phase boundary 
between the collapsed and the spiral phase is still unknown. As mentioned in Sec. [II B| , 
collapsed polymers are different from dense polymers, for which C = 4. Secondly, the 
crossover exponent \1' introduced in Eq. @ needs to be determined. Only if the ^-point 
is a special point in the phase diagram, the number of parallel contacts increases with A^* 
at the ^-point. Otherwise, it increases logarithmically in N, i.e., \E' vanishes. Finally, 
the transition from the swollen phase to the spiral phase is poorly understood. Although 
numerical results indicate a first-order transition we cannot rule out a continuous 

transition. The assumption in that the spiral has no entropy cannot be upheld, as 
we have seen in our study of the globule-to-spiral transition. It is quite possible that the 
swollen-to-spiral transition has some analogy with the adsorption of a polymer at a wall, 
which is known to be continuous. The "wall" in our case is the surface of a spiral arm, and 
the thickness of the adsorbed polymer layer corresponds to the thickness of a spiral arm. 
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FIG. 1. An oriented self-avoiding path on a square lattice 
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FIG. 2. Qualitative phase diagram of interacting self-avoiding walks 
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FIG. 3. An interacting self-avoiding walk on the Manhattan lattice at the collapse transition 





FIG. 5. Winding angle distribution at the collapse transition on the Manhattan lattice 
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FIG. 6. A spiral composed of globules 
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FIG. 7. Value of the constant C in most parts of the phase 
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